from subprocess import call

f = open("LongRead_DIV", 'w')

for i in range(12,21):
	print "\nNow simulating seed length: %i" % i;
	f.write('\n***SeedLength: %i\n' % i);
	f.flush();
	call (["./mapSimulator", "human_g1k_v37.fasta", ("human_g1k_v37.fasta.ht." + str(i) ), '-l', "testRead"], stdout=f);

